Ab initio calculations of phonon spectra in ATiOs perovskite crystals 
(A = Ca, Sr, Ba, Ra, Cd, Zn, Mg, Ge, Sn, Pb) 
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The phonon spectra of calcium, strontium, barium, radium, cadmium, zinc, magnesium, germa- 
nium, tin, and lead titanates with the perovskite structure are calculated from first principles within 
the density functional theory. By analyzing the unstable modes in the phonon spectra, the possible 
lattice distortions are determined and the energies of the corresponding phases are calculated. From 
analyzing the phonon spectra, force constants, and eigenvectors of TO phonons, a conclusion is 
drawn on the origin of the ferroelectricity in considered crystals. It is shown that the main fac- 
tors determining the possible off-centering of atoms in the A position are the geometric size and 
electronic configuration of these atoms. 
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I. INTRODUCTION 

Crystals of the perovskite family are well-known ma- 
terials undergoing various structural distortions with de- 
creasing temperature. When the character of these dis- 
tortions is ferroelectric, a number of physical properties 
of these crystals (dielectric constant, piezoelectric coeffi- 
cients, etc.) becomes anomalously large. For this reason, 
these materials have found wide application in modern 
electronics. 

The problem of further optimization of the ferroelec- 
tric properties requires a deeper understanding of the mi- 
croscopic mechanisms responsible for the appearance of 
ferroelectricity and the ferroelectric properties. In solv- 
ing this problem, very useful information can be obtained 
from the ab initio calculations which have already made 
a significant contribution to the understanding of the fer- 
roelectric phenomena in the perovskite crystalsj^^— 

In discussing the properties of ferroelectrics, it is very 
important to understand whether these properties result 
from a collective displacement of atoms in a lattice (dis- 
placive phase transition) or they are due to specific fea- 
tures of certain constituent atoms in a crystal (order- 
disorder phase transition). This problem arises, in par- 
ticular, in discussing the nature of the phase transitions 
that occur in incipient ferroelectrics doped with certain 
impurities. 

Earlier studies of titanates with the perovskite 
structure have dealt mainly with four compounds 
(CaTiOa^iiiHiii SrTiOg^iii^ BaTiO^i^-^^ and 
PbTiOa^"— i^i^) and their solid solutions. However, a 
comparison of the results obtained in these studies 
is hampered by different calculation techniques used 
in these investigations and different methods used to 
construct atomic (pseudo)potentials. 

The aim of this work is to carry out first-principles 
calculations of the phonon spectra of ten ^TiOa crystals 
with the perovskite structure and determine the struc- 
ture of the most stable phases. In order to test our ap- 
proach, we first apply it for calculating the properties 



of the four above-mentioned systems. Then we predict 
the properties of poorly investigated or hypothetical per- 
ovskite compounds RaTiOa, CdTiOa, MgTiOs, ZnTiOa, 
SnTiOs, and GeTiOs.^^ From comparison of the results 
obtained in a unified way for a large number of related 
materials, conclusions are made about the relation of the 
structural distortions in ^TiOs crystals to the size and 
electronic structure of the A atom. From the analysis 
of the on-site force-constant matrix and the TO-phonon 
eigenvectors, we draw conclusions about the origin of the 
ferroelectricity in these materials and find conditions un- 
der which the ferroelectric phenomena can be associated 
with off-center atoms. 



II. CALCULATION TECHNIQUE 

The calculations were performed using the ABINIT 
software^ - based on the density functional theory, pseu- 
dopotentials, and the plane-wave expansion of the wave 
functions. The exchange-correlation interaction was de- 
scribed in the local density approximation (LDA)jii The 
used pseudopotentials were optimized separable nonlo- 
cal pseudopotentials,= which were constructed using the 
OPIUM code and to which the local potential correc- 
tion^^ was added in order to improve their transferability. 
For elements with atomic numbers Z < 46, the construc- 
tion was performed in a non-relativistic way; for other 
elements, a scalar-relativistic approximation was used. 
Table |T] lists the parameters used for constructing pseu- 
dopotentials. The local potential was the s potential ex- 
cept for the oxygen atom, for which the local d potential 
was used. The parameters of pseudopotentials were finely 
adjusted by comparing the calculated and experimental 
values of the lattice parameters for a number of oxides 
and sulfides of the elements. 

The lattice parameters and the equilibrium atomic po- 
sitions in the unit cell were obtained by minimizing the 
Hellmann-Feynman forces acting on the atoms (<10~^ 
Ha/Bohr), with the total energy being calculated self- 
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consistently with an accuracy of better than 10~^° Hai^S 
In the calculations, particular attention was paid to the 
convergence of the results with respect to the kinetic en- 
ergy cut-off for plane waves and the density of the fc-point 
mesh used in the integration over the Brillouin zone. For 
all the calculated properties presented below, the conver- 
gence was attained at an energy cut-off of 30 Ha and the 
8x8x8 fc-point mesh constructed according to Ref . [2l|. 

The Born effective charges Z*, optical dielectric con- 
stant Eoo, elastic moduli C^j, bulk modulus B, force- 
constant matrix $.y , and phonon spectra were calculated 
using the density-functional perturbation theoryi22ri2^ 
The phonon frequencies were calculated exactly at five 
points of the Brillouin zone (F, X, M, R, and the A 
point located halfway between the F and R points), and 
then the phonon spectrum was computed over the entire 
Brillouin zone using the interpolation technique i^^i^^ 



III. TESTING OF THE CALCULATION 
TECHNIQUE 

The correctness of the described approach was tested 
by comparing the calculated lattice parameters, sponta- 
neous polarizations, and phonon spectra with available 
experimental data and calculations performed by other 
authors for well-studied CaTiOs, SrTiOs, BaTiOs, and 
PbTiOa compounds. 

The lattice parameters corresponding to a minimum 
total energy of the crystals are given in Table [TTl The 
obtained values agree well with the experimental data^^ 
if one takes into account that the LDA usually slightly 
underestimates the lattice parameter. An analysis of 
the relative energies (per formula unit) of low-symmetry 
phases (Table IIII[) shows that for barium titanate, the 
most energetically favorable phase is the rhombohedral 
one and for lead titanate, it is the tetragonal phase. 
The calculated values of the c/a ratio for tetragonal 
BaTiOa and PbTiOa are close to the experimental val- 
ues (Table HH . For CaTiOa, the most energetically fa- 
vorable phase is the orthorhombic Pbnm phase and for 
SrTiOs, it is the tetragonal I A/ mem phase. The values of 
spontaneous polarization calculated by the Berry phase 
method^ are 0.26, 0.31, and 0.89 C/m^ for tetragonal 
BaTiOa, rhombohedral BaTiOs, and tetragonal PbTiOs, 
respectively. These values are close to the experimental 
data (0.26, 0.33, 0.75 C/m^, Ref.IIg). 

The calculated optical-phonon frequencies are also in 
good agreement with available experimental data and 
calculations performed by other authors (Tables IIV| W\ . 
The imaginary frequencies in Tables llVl andlVlcorresDond 
to unstable modes, the squares of whose frequencies are 
negative. 

The good agreement of our calculations with the ex- 
perimental data and the calculations of other authors for 
CaTiOg, SrTiOs, BaTiOg, and PbTiOg suggests that the 
proposed approach can be used to predict the properties 
of poorly studied and hypothetical2i titanates with the 



perovskite structure and to discuss the factors causing 
the appearance of ferroelectricity in them. The proper- 
ties of these crystals calculated for the theoretical lattice 
parameter (i.e., the value corresponding to a minimum of 
the total energy) are given in Tables IIIII to IVIIl Table IIIII 
presents the energies of different low-symmetry phases 
measured relative to the energy of the parent cubic phase. 
Table HVl gives the frequencies of optical phonons (three 
infrared-active F15 modes and one infrared-inactive F25 
mode) in the cubic phase. Table |V] gives the lowest 
phonon frequencies at high-symmetry points of the Bril- 
louin zone. The values of the Born effective charges and 
optical dielectric constants for the cubic phase of ^TiOa 
compounds are given in Table IVIl Finally, the elastic 
moduli of several crystals are presented in Table IVIIl 

The phonon dispersion curves along some high- 
symmetry directions calculated for the cubic phase of ten 
titanates are shown in Fig.[TJ The imaginary phonon fre- 
quencies associated with the structural instability of the 
crystals are represented in Fig. [T] by negative values. 



IV. RESULTS 

As follows from Fig. [1] unstable optical modes of dif- 
ferent symmetry are present in the phonon spectra of all 
studied titanates. At first, we discuss the phonon spectra 
of well-studied compounds. 

A specific feature of the dispersion curves of SrTiOa 
is that among three unstable phonons at the F, i?, and 
M points, the most unstable phonon is that at the R 
point (the R25 mode)^ and that the phonon frequency 
depends only weakly on the wave vector along the R-M 
line. As shown in Ref. d, in the real space, these unstable 
phonons with the wave vectors near the edges of the cubic 
Brillouin zone correspond to the rotation of the oxygen 
octahedra linked together by shared vertices, with the 
correlation length for rotations being three to five lattice 
periods. Thus, the R25 mode and the M3 mode (which 
corresponds to the unstable phonon at the M point) de- 
scribe the structural distortions associated with the ro- 
tation of the octahedra. The unstable F15 phonon mode 
at the F point corresponds to the ferroelectric distortion 
of the structure. 

The 7?25 mode is triply degenerate, and the distor- 
tions described by the three-component order parameters 
(77, 0, 0), (77, 77, 0), and (77, -q, rf) lead to low-symmetry 
phases with the space groups lA/mcm, Imma, and i?3c, 
respectively. The M3 mode is nondegenerate, and its 
condensation lowers the crystal symmetry to PA/mbm. 
The triply degenerate Fis mode results in the well-known 
space groups PAmm, Amm2, and R3m. From compar- 
ing the energies of these phases (Table IIIip . it follows 
that the lowest-energy phase in SrTiOa is the I A/ mem 
one, into which strontium titanate transforms with de- 
creasing temperature. The instability of the ferroelectric 
mode in this compound is not sufficiently strong for fer- 
roelectricity to occur. 
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TABLE I. Electronic configurations of atoms and parameters used to construct pseudopotentials: r^, rp, and are tlie core 
radii of the pseudopotentials for the s, p, and d projections; qs, Qp, and qd are the cut-off wave vectors used to optimize 
pseudopotentials; and rmin, rmax, and VIoc are the limits and the depth of the correcting local potential (parameter values are 
in atomic units, and energy is given in Ry). 



Atom 


Configuration 




rp 




qs 


<lv 


<ld 




^max 




Ca 


3s23p«3d"4s'' 


1.46 


1.68 


1.82 


7.07 


7.07 


7.27 


0.01 


1.40 


1.6 


Sr 




1.68 


1.74 


1.68 


7.07 


7.07 


7.07 


0.01 


1.52 


1.5 


Ba 




1.85 


1.78 


1.83 


7.07 


7.07 


7.07 


0.01 


1.68 


1.95 


Ra 




1.84 


1.73 


1.98 


7.8 


7.8 


7.8 


0.01 


1.68 


-1.3 


Mg 


2s22p«3s03pO 


1.50 


1.88 




6.7 


8.1 




0.01 


1.0 


-0.84 


Zn 


3di»4s°V 


1.82 


1.82 


2.00 


7.07 


7.07 


7.47 


0.01 


1.60 


2.5 


Cd 


4di°5s°5p° 


2.04 


2.18 


2.10 


7.07 


7.07 


7.07 





1.88 


-1.6 


Ge 


3d"4si-54p°-5 


1.68 


1.68 


1.96 


7.07 


6.0 


7.77 


0.01 


1.58 


0.48 


Sn 


4di«5s2 5pO 


2.14 


2.08 


2.18 


7.07 


7.07 


7.07 


0.01 


1.90 


0.64 


Pb 


Sd^^Gs^epO 


1.72 


1.98 


1.80 


6.05 


5.52 


7.17 


0.1 


1.43 


1.6 


Ti 


3s23p83d04s0 


1.48 


1.72 


1.84 


7.07 


7.07 


7.07 


0.01 


1.41 


2.65 







1.40 


1.55 


1.40 


7.07 


7.57 


7.07 









TABLE II. Comparison of calculated and experimental lattice parameters of different phases of ATiOa compounds (experimental 
data are obtained at 300 K, unless otherwise specified). 



Compound 


Space group 


Source 


Lattice parameters 


CaTiOs 


Pbnm 


This work 


a = 5.3108, b = 5.4459, c = 7.5718 A 






Exp.^ 


a = 5.3670, b = 5.4439, c = 7.6438 A 


SrTiOa 


Pm3m 


This work 


a = 3.8898 A 






Exp.^ 


a = 3.905 A 




1 4/ mem 


This work 


a = b = 5.4680, c = 7.8338 A 






Exp.^ 


a^b = 5.510, c = 7.798 A (20 K) 


BaTiOa 


Pm3m 


This work 


a = 3.9721 A 






Exp.^ 


a = 3.996 A (393 K) 




PAmm 


This work 


a = 3.9650, c = 4.0070 A, c/a = 1.0106 






Exp.^ 


a = 3.9920, c = 4.0361 A (293 K) 




Amm2 


This work 


a = 3.9620, b = 5.6384, c = 5.6484 A 






Exp.^ 


a = 3.990, b = 5.669, c = 5.682 A (263 K) 




R3m 


This work 


a = 3.9817 A, Q = 89.933° 






Exp.^ 


a = 4.001 A, Q = 89.85° (105 K) 


PbTiOs 


PAmm 


This work 


a = 3.8858, c = 4.1151 A, c/a = 1.0590 






Exp.2^ 


a = 3.904, c = 4.152 A 



In CaTiOa, in addition to the instabilities found above, 
three weak antiferroelectric-type instabilities associated 
with the X^, X^, and modes arise and the R- 
M segment of the phonon spectrum becomes practi- 
cally dispersionless (compare mode energies in Table IVT) . 
In the latter case, as shown by theoretical calcula- 
tions, ^4 the simultaneous condensation of two unstable 
i?25 and M3 modes results in the formation of the Pbnm 
phase having the lowest energy among the possible dis- 
torted phases ( Table HIH .'^'^ The transition from the high- 
temperature PmSm phase to the low-temperature Pbnm 
phase can occur through one of three intermediate phases 
{PA/mbm, lA/mcm, and -R3c), whose energies are 0.11- 



0.17 eV higher than the energy of the Pbnm phase.— 
The ferroelectric PAmm and i?3m phases in CaTiOs have 
much higher energies and never arise. As for the weakly 
unstable X'^, X^, and Mg modes, they are doubly de- 
generate and the distortions described by the order pa- 
rameters (ry, 0) and (r/, 77) result in the Pmma, Cmcm, 
and Cmmm phases. The energy gain resulting from the 
transformation into these phases, however, does not ex- 
ceed 7 meV. 

The phonon spectrum of BaTiOa significantly differs 
from the spectra discussed above by the absence of in- 
stability at the R point and the appearance of highly 
unstable and M3 modes (at the X and M points. 
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TABLE III. Relative energies of different low-symmetry phases of ATiOs compounds (energies of the most stable phases are 
in boldface). 



Compound 


Unstable 


Space 


Energy 


Unstable 


Space 


Energy 


Unstable 


Space 


Energy 




mode 


group 


(meV) 


mode 


group 


(meV) 


mode 


group 


(meV) 


MgTiOs 


r25 


Pim2 


-125" 


Mi 


Cmmm 


-417" 


TlSi r25 


Amm2 


-1380" 




X3 


P4:2/mmc 


-147" 


X'5 


Pmma 


-500 


R2S + Ms 


Cmcm 


-1658" 






M/mmm 


-228" 


T25 


R32 


-686" 


R25 


R3c 


-1727 






RSm 


-289" 


Tib 


RSm 


-695 




Imma 


-1764" 






Cmcm 


-304 


Tis 


PAmm 


-1028 


R25 + Ms 


Pbnm 


-1992 




Mi 


Pmma 


-344" 


Ms 


PA/mbm 


-1107 










Ml, 


P'i/nmm 


-417" 


R25 


I A/ mem 


-1111 








CaTiOa 




Cmcm 


-0.0" 


Mi 


Cmmm 


-6.7 


R25 


lA/mcm 


-365 




X5 


Pmma 


-0.0" 


Tis 


R3m 


-73.7 


R25 


R3c 


-385 






Cmcm 


-0.6 


Tib 


Amm2 


-85.4" 


R25 + Ms 


Cmcm 


-404" 




xk 


Pmma 


-0.9 


Tib 


PAmm 


-123 


R26 


Imma 


-412" 




Mi, 


Pmma 


-5.0 


Ms 


PA/mbm 


-321 


R25 + Ms 


Pbnm 


-497 


SrTiOs 


ri5 


PAmm 


-0.71 


Ms 


PA/mbm 


-9.45 


R26 


I A/ mem 


-30.9 






Amm2 


-0.75" 


R25 


R3c 


-27.5 












R3m 


-0.75 


R25 


Imma^ 


-28.9 








BaTiOg 


M3 


PA/nmm 


-0.31 


X5 


Cmcm 


-1.45 


Tib 


Amm2 


-7.4 






Pmma 


-1.22 


Tis 


PAmm 


-5.6 


Tib 


R3m 


8.1 


RaTiOs 


M3 


PA/nmm 


-11.1 


X5 


Cmcm 


-16.9 


Tib 


Amm2 


-28.5 




Xb 


Pmma 


-14.2 


Tib 


PAmm 


-21.8 


Tib 


R3m 


-29.7 


CdTiOs 


R16 


lA/mmm 


-24" 


Xi 


Pmma 


-160 


R25 


I A/ mem 


-912 




Rir, 


R3m 


-30" 


Tis 


R3m 


-245 


Ms 


PA/mbm 


-920 




Xs 


P4:2/mmc 


-45 


X5 


Cmcm 


-282 


R25 + Ms 


Cmcm 


-1157" 




X5 


Cmcm 


-66" 


Tib 


PAmm 


-340 


R25 


R3c 


-1197 




xk 


Pmma 


-104" 




Amm2 


-412 


-/?15, R2S 


Imma 


-1202" 






PAm2 


-134 


r2B 


R32 


-486 


R25 + Ms 


Pbnm 


-1283 


ZnTiOs 


X3 


PA^/mmc 


-171" 


Mi 


Cmmm 


-564" 


R26 


I A/ mem 


-1443 






Plm2 


-341 


M2 


PA/nmm 


-688" 


Ms 


PA/mbm 


-1449 




fil5 


M/mmm 


-375" 


xi 


Pmma 


-752" 


r25 


R32 


-1486 






Pmma 


-447 


Xf. 


Cmcm 


-867 


R2S + Ms 


Cmcm 


-2036" 




^5 


Cmcm 


-449" 


Tib 


R3m 


-868 


R15, R25 


Imma 


-2215" 






RSm 


-465" 


Tis 


PAmm 


-1104 


R25 


R3c 


-2271 




ML, 


Pmma 


-555" 


TlS, r2B 


Amm2 


-1254 


i?25 + Ms 


Pbnm 


-2312 


GeTiOs 




Plm2 


-0.5" 


Ms 


PA/mbm 


-AAA 


R20 + Ms 


Cmcm 


-733" 




r25 


RA2 


-1.5" 


R25 


I A/ mem 


-455 


R25 + Ms 


Pbnm 


-810 




Rl5 


1 4 1 mmm 


-314" 


R15 


R3m 


-501" 


Tib 


PAmm 


-854 




X', 


Pmma 


-328 


i?25 


B3c 


-589 


riv r2r, 


Amm.2 


-881" 




x^ 


( 'mciii 


-12« 


Rl:,, R2:, 


I ninia 


-70i" 


I'lr, 


R'im 


-1053 


SnTiOa 


Rl5 


I A/ mmm 


-1.1" 


Ms 


PA/mbm 


-57 


R15, R25 


Imma^ 


-84 




R16 


i?3m 


-1.7" 


R25 


lA/mcm 


-67 


Tib 


R3m 


-240 




x^ 


Pmma 


-21 


R25 


R3c 


-74 


Tib 


Amm2 


-242" 




x^ 


Cmcm 


-23 


R25 + Ms 


Cmcm 


-78" 


Tib 


PAmm 


-291 


PbTiOs 


Ms 


PA/mbm 


-10.1 


R25 


R3c 


-21.6 


ri5 


Amm2 


-70.4" 




R2-. 


I A/ mem 


-19.6 


R25 


Imma^ 


-22.2 


Tib 


PAmm 


-84.4 




R2r, -T- Ms 


Cmcm 


-19.7" 


Tib 


R3m 


-66.3 









"Results of additional calculations. 

''The Pbnm structure found in earlier calculations relaxes to the Imma one. 
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TABLE IV. Frequencies of optical phonons at the F point in 
the cubic phase of ylTiOs compounds (in cm^^). 



Compound 


Source 


1 Vjl 


1 V^Z 


T03 LOl 






i 25 


MgTiOa 


this work 


260^ 


151 


649 


106j 


372 


905 


I9I2 


UaiiUa 


this work 




lib 


607 


122 


4U ( 


OO ( 


yo 




theor.— 


1532 


188 


610 


133 


427 


866 






theor^^ 


1 Ar\A 


zUU 


625 


136 




od4 


ioU 


brilUs 


this work 


Do^ 


IDZ 


549 


152 


A OO 

4zo 


'700 


ono 
2U2 




theor.— 


4i2 


loo 


546 


158 


A A 

4o4 








tneor.— 


1 nn^' 


1 K 1 

lot 


522 


146 


A on 

4oy 


(^Ol 


01 n 








i ( 


544 


172 


A 

4 ^ 


/yo 




rsaiiUa 


this work 


iOi2 


i ( 


471 


172 


4oy 


DE5o 


zoy 








1 77 


468 


173 




( OO 






tneor." 




IDD 


455 


162 


4o4 


DO ( 


ZDD 








1 ft! 
ioi 


487 


180 


4Do 


71 7 


oUD 


riailvJ3 


this work 




1 79 


444 


166 


A 'iA 


Ooo 




ZnTiOa 


this work 


2402 


76 


645 


105i 


316 


815 


3532 


CdTiOg 


this work 


187z 


97 


616 


34 


353 


820 


2312 


GeTiOg 


this work 


2472 


122 


583 


68i 


356 


762 


49i 


SnTiOs 


this work 


185z 


126 


505 


80 


375 


689 


183 


PbTiOa 


this work 


1502 


116 


499 


96 


394 


693 


202 




theor « 


1442 


121 


497 


104 


410 


673 






theor.— 


1822 


63 


447 


47 


418 


610 





Data for the tetragonal phase. 



respectively) corresponding to the antiferroelectric dis- 
tortions of the structure into the Pmma, Cmcm, and 
PA/nmm phases. Among the unstable modes, the fer- 
roelectric Fis mode is the most unstable and it deter- 
mines the distortions observed in the crystal (the ener- 
gies of three above-mentioned antiferroelectric phases are 
higher than that of the polar PAmm phase). We note 
that the weaker phonon instability in our calculations as 
compared to the results of Ref. |^ is because our calcula- 
tions were performed for the theoretical lattice parame- 
ter, whereas the calculations of Ref. were performed for 
its experimental value. To illustrate the influence of the 
lattice strain on the phonon spectrum, the low-energy 
part of this spectrum calculated for the same lattice pa- 
rameter as in Ref. |^ is shown in Fig. [T] by the dashed 
line. 

The weak dependence of the unstable TO phonon fre- 
quency on the wave vector along the F-X-M-F path for 
vibrations polarized along the fourfold axes of the cubic 
structure was first discovered in Ref. J5. This dependence 
gives evidence for the dominant role of the vibrations in 
the ...-0-Ti-O-... linear chains oriented along these axes 
and for the weak interaction between adjacent chains. 

A comparison of the phonon spectra of BaTiOs and 
RaTiOa shows that these spectra are very similar. In ra- 
dium titanate, the most unstable mode is the F15 mode 
whose frequency depends only weakly on the wave vector 
along the F-AT-M-F path. The instability of phonons in 



TABLE V. Lowest phonon frequencies at high-symmetry 
points of the Brillouin zone in the cubic phase of j4Ti03 com- 
pounds (in cm~^). 



Compound 


Source 


r 


X 


M R 


A 


MeTiO'^ 

ivj.g^ _i_ J. V y ^ 


this work 


260j 


190i 


314i 315j 


246j 


CaTiOa 


this work 


165i 


32j 


215i 226j 


122j 




theor.— 


140i 


20 


207i 219j 




SrTiOa 


this work 


68i 


98 


86j 119i 


100 


BaTiOa 


this work 


151i 


96i 


59i 134 


105 




theor.** 


219i 


189i 


167i 128 




RaTiOs 


this work 


212i 


182i 


158i 110 


87 


ZnTiOs 


this work 


353i ° 


319i 


437i 424i 


337i 


CdTiOa 


this work 


2311" 


184i 


333i 328i 


265i 


GeTiOs 


this work 


247i 


148i 


254i 251j 


201j 


SnTiOa 


this work 


185i 


56i 


144i 148j 


97i 


PbTiOs 


this work 


150i 


30 


96i 113* 


15i 




theor.— 


182i 


31i 


35i 145i 


58i 


"The mode. 










TABLE VL Effective charges and optical dielectric constant 


for the cubic phase of ATiOs compounds. 






Compound 


zi 




Zoi_ 


■^011 


f 00 


MgTiOs 


2.537 


7.773 


-2.026 


-6.258 


7.01 


CaTiOa 


2.579 


7.692 


-2.085 


-6.101 


6.84 


SrTiOs 


2.561 


7.725 


-2.099 


-6.088 


6.87 


BaTiOa 


2.738 


7.761 


-2.186 


-6.128 


7.28 


RaTiOa 


2.764 


7.789 


-2.181 


-6.192 


7.42 


ZnTiOs 


3.233 


8.257 


-2.427 


-6.637 


11.64 


CdTiOa 


3.040 


8.052 


-2.300 


-6.493 


9.35 


GeTiOs 


4.460 


7.572 


-2.860 


-6.314 


10.49 


SnTiOs 


4.255 


7.529 


-2.745 


-6.294 


10.18 


PbTiOa 


3.931 


7.623 


-2.635 


-6.283 


9.34 



TABLE VII. Elastic moduli of the cubic phase of ^TiOa com- 
pounds (in GPa). 



Compound 


Source 


Cii 


C12 


C44 


B 


CaTiOa 


this work 


388 


100 


91 


196 




theor.* 


407 


96 


101 


200 


SrTiOs 


this work 


373 


103 


108 


193 




theor.— 


388 


104 


117 


199 




27 

exp.— 


316-348 


101-103 


119-124 


174-183 


BaTiOa 


this work 


338 


110 


123 


186 




theor.— 


329 


117 


130 


188 


RaTiOs 


this work 


319 


112 


126 


181 


PbTiOs 


this work 


336 


127 


95 


197 




theor.— 


334 


145 


100 


208 




theor.— 


320 


141 


375? 


203 
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RaTiOa is even more pronounced than in BaTiOa. Tak- 
ing into account the calculated energies of the distorted 
phases (Table IIII|) , we can suppose that RaTiOs is also 
a ferroelectric and that, as the temperature decreases, it 
undergoes three successive phase transitions as barium 
titanate does. The temperatures of these transitions are 
likely to exceed those in barium titanate. The values 
of the spontaneous polarization calculated by the Berry 
phase method in RaTiOs are also higher than those in 
BaTiOa; they are 0.36 G/vo? in the tetragonal phase and 
0.41 C/w? in the rhombohedral phase. The calculated 
elastic moduli of the cubic RaTiOs are given in Table IVlIl 

Now, we consider the phonon spectrum of CdTiOs. 
This spectrum has a number of unstable modes at the 
X, Af, and R points (Xg, Xg, X'^, M3, M^, M^, i?25, 
i?i5 modes) and two unstable modes at the T point. It is 
surprising that the main instability at the V point is due 
to the mode associated with the deformation of the 
oxygen octahedron rather than to the ferroelectric Fig 
mode (see mode frequencies in Table ITVT) . This deforma- 
tion can result in the formation of the P4m2, Ammfl, 
and i?32 phases, depending on the number of nonzero 
components of the order parameter. '^^ The energy of the 
most stable of these phases (i?32, Table UlI)) is lower than 
the energy of the polar phases. Because of the qualita- 
tive similarity between the phonon spectra of calcium 
and cadmium titanates and between the eigenvectors of 
their unstable modes at the R and M points, CdTiOa 
can be considered as an analog of CaTiOs characterized 
by an even higher instability. Therefore, at room tem- 
perature, the structure of its nonpolar phase is Pbnm, as 
in the case of CaTiOg. This conclusion was also made 
in Ref. According to our data, the energy of this 
phase is lower than that of the cubic phase by 1.28 eV 
(Table IIII[) , which is somewhat higher than the value ob- 
tained in Ref. 36 (0.8 eV) and Ref. M (0.91 eV). 

Although the ferroelectric instability associated with 
the mode is not very important in the cubic CdTiOs, 
it is known that this instability exists in the Pbnm phase 
and results in the ferroelectric phase transition at 80 K. 
The first-principles calculations of the properties of the 
orthorhombic Pbnm phas^^ did not found a stable fer- 
roelectric distortion in it. In contrast, our calculations of 
the phonon spectrum at the T point of the orthorhombic 
CdTiOa reveal two unstable i3i„ and modes resulting 
in the formation of polar Pb2im and Pbn2i phases, re- 
spectively. These lattice distortions have been observed 
in X-ray studies of cadmium titanate at low tempera- 
tures i^i^ The properties of these phases will be consid- 
ered in a separate paper. 

The phonon spectrum of ZnTiOg is qualitatively simi- 
lar to that of CdTiOs, but it has an additional weak un- 
stable mode and is even less stable. The mode in 
it is also less stable than the ferroelectric Fis mode (Ta- 
ble lIV[) . However, as the most unstable modes in ZnTiOs 
are the R25 and M3 ones, the Pbnm phase has the lowest 
energy (Table HIX)) . The calculations of the phonon spec- 
trum at the r point of the orthorhombic zinc titanate re- 



veal two unstable Biu and modes, which can cause 
the formation of the same polar phases as in cadmium 
titanate. 

The phonon spectrum of MgTiOs is intermediate be- 
tween those of zinc titanate and calcium titanate. It 
also has unstable Fis and modes, but in magnesium 
titanate the ferroelectric Fig mode is less stable (Ta- 
ble |IV|. Nevertheless, as the phonons at the R and M 
points have the lowest frequency, the Pbnm phase is the 
most energetically favorable (Table HIT)) . The calculations 
of the phonon spectrum at the F point in the orthorhom- 
bic magnesium titanate reveal one unstable Biu mode, 
which can cause the Pbnm — >■ Pbn2i ferroelectric phase 
transition. 

Finally, we discuss the phonon spectra of PbTiOa, 
SnTiOs, and GeTiOa. The ferroelectric instability in 
these three compounds is associated with the F15 mode 
which competes with the unstable R25 and M3 modes. 
From comparing the energies of different distorted phases 
(Table |TTT|, it follows that even in GeTiOs, in which the 
unstable phonons at the F, R, and M points are close 
in frequency, the ferroelectric instability is dominant. 
The calculated spontaneous polarization is 1.28 C/m^ 
in the tetragonal SnTiOa and 1.37 C/m^ in the rhombo- 
hedral GeTiOg. Our value of the spontaneous polariza- 
tion in SnTiOg is significantly higher than the value of 
0.73 C/m^ obtained in Ref. 40:. Among the perovskite 
compounds studied to date, GeTiOs is likely to have the 
highest spontaneous polarization. 

It should be noted that in PbTiOs and SnTiOs, 
the most energetically favorable phase is the tetragonal 
PAmm phase, whereas in GeTiOg it is the rhombohedral 
R3m phase. At the same time, the lattice strain in the 
tetragonal GeTiOa (c/a — 1.1821) is much higher than in 
lead titanate (c/a = 1.0590). This result calls into ques- 
tion the conclusion^ that the stabilization of the tetrago- 
nal phase in PbTiOg is due to a high lattice strain (large 
c/a ratio). 



V. DISCUSSION 

As follows from Fig.[Tl the phonon spectra of all ATiO^ 
perovskite crystals studied in this work are characterized 
by several unstable modes, one of which is always the fer- 
roelectric F15 mode. In the case where the R25 and M3 
modes competing with it have a lower frequency, the crys- 
tal undergoes distortions — the octahedra rotations, and 
its symmetry is lowered to I A/ mem or Pbnm. The ten- 
dency toward such structural phase transitions increases 
with decreasing the A atom size. 

From analyzing the characteristics of the F15 mode, 
one can draw a conclusion on the nature of the ferroelec- 
tric instability in studied crystals. As mentioned above, 
the dispersion law of this mode in BaTiOa and RaTiOs 
indicates strongly correlated motion of atoms along the 
...-0-Ti-O-... chains. The analysis of the F15 phonon 
eigenvectors (Table IVIIip shows that the A atoms do not 
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TABLE VIII. Eigenvectors of the dynamic matrix for an un- 
stable TOl phonon at the F point in the cubic phase of ylTiOs 
compounds. 



Compound 








a^oii 


MgTiOa 


+0.6828 


+0.1831 


-0.4800 


-0.1985 


CaTiOg 


+0.5693 


+0.2391 


-0.5225 


-0.2696 


SrTiOa 


+0.3434 


+0.3852 


-0.5372 


-0.3956 


BaTiOa 


+0.0299 


+0.6734 


-0.3561 


-0.5404 


RaTiOa 


+0.0051 


+0.6750 


-0.2841 


-0.6188 


ZnTiOs 


+0.5167 


+0.1889 


-0.5655 


-0.2403 


CdTiOa 


+0.4012 


+0.2358 


-0.5919 


-0.2875 


GeTiOs 


+0.5367 


+0.1382 


-0.5573 


-0.2677 


SnTiOa 


+0.4177 


+0.2123 


-0.5670 


-0.3709 


PbTiOs 


+0.2973 


+0.2865 


-0.5675 


-0.4305 


TABLE IX. 


Diagonal elements of the 


on-site force-constant 


matrix 


'0, 0) for A and Ti atoms 


in the cubic phase of 


^TiOa compounds (in Ha/Bohr^). 






Compound 




A atom 




Ti atom 


MgTiOa 




-0.0109 




+0.1431 


CaTiOs 




+0.0163 




+0.1370 


SrTiOs 




+0.0445 




+0.1196 


BaTiOa 




+0.0755 




+0.0873 


RaTiOa 




+0.0856 




+0.0750 


ZnTiOa 




-0.0229 




+0.1072 


CdTiOa 




-0.0008 




+0.1113 


GeTiOa 




-0.0150 




+0.0949 


SnTiOa 




+0.0132 




+0.0786 


PbTiOa 




+0.0269 




+0.0803 



contribute niuch to the motion which is mainly deter- 
mined by the out-of-phase Ti and 0|| displacements. As 
the size of the A atoms decreases, their contribution to 
the motion increases and becomes dominant, whereas the 
contribution of the Ti atoms decreases and the out-of- 
phase motion involves now not 0|| but displacements. 
Thus, in crystals with small A atoms, the ferroelectric 
mode is determined by the out-of-phase displacements of 
the A atom and the cuboctahedron of oxygen atoms. 

The values of the diagonal elements of the on-site force- 
constant matrix $xx(0, 0) for A and Ti atoms are pre- 
sented in Table IIXI These matrices are defined in terms 
of the restoring force acting on an atom displaced from its 
site, with the other atoms fixed at their lattice sites. To 
determine the on-site force constants from the force con- 
stants calculated using the ABINIT software for a sublat- 
tice displaced as a whole, the force constants calculated 
on a regular mesh of wave vectors were averaged-^^^iSli^ 
Positive values of the on-site force constants indicate that 
the position of an atom at its site is stable, while the 
negative values indicate that the atom is off-center. As 



0.10 

0.08 

^ 0.06 
o 
ffl 

0.04 
S 0.02 









Ra 








Ba 


- 






sr 








X jm Pb 




Mg 

Zn 

1 




1 1 


0.6 


0.8 


1.0 


1.2 1.4 



FIG. 2. Dependence of the $12, (0,0) matrix element for the 
A atom on its ionic radius. 



follows from Table 11X1 the off-centering of the A atoms 
should be observed in ^TiOs perovskites for A = Mg, 
Zn, Cd, and Ge. The Sn, Ca, and Pb atoms are fairly 
close to the stability limit against the transition to an 
off-center position. 

It should be recalled that in this work, the calculations 
were performed for the theoretical lattice parameter (cor- 
responding to a minimum of the total energy). As the 
systematically underestimated lattice parameters in the 
LDA weakens the ferroelectric instability, many authors 
perform calculations for the experimental lattice parame- 
ters. In order to estimate the influence of this systematic 
error, we carried out a computer simulation which showed 
that an increase in the lattice parameter by 1% (which 
is a typical LDA error) decreases $a:£c(0, 0) in PbTiOs by 
0.006 Ha/Bohr2 for the A atom and by 0.016 Ha/Bohr^ 
for the Ti atom. As a result, the atoms with $a;a;(0,0) 
close to the stability limit against the transition to an 
off-center position can actually be off-center. Perhaps, 
this occurs in lead titanate, as indicated by extended X- 
ray-absorption fine-structure (EXAFS) studies^i 

The main parameter determining the tendency of the 
A atom in ^TiOs compounds to become off-center is its 
atomic size. However, if one plots the dependence of 
*i'a;a(0, 0) on the ionic radius of the A atom (Fig. [5]), it 
becomes clear that the curve for Zn and Cd atoms, as well 
as the curve for Ge, Sn, and Pb atoms, are different from 
that for the "main" Mg-Ca-Sr-Ba-Ra series. This differ- 
ence is likely a consequence of different configuration of 
the outer electron shell for the atoms. This configuration 
is rfio for 2n and Cd; d^^ s"^ for Ge, Sn, and Pb; and sV 
for the atoms of the main series. The difference in the 
properties of these groups of atoms is clearly manifested 
in the Born effective charges of the A atoms (Table IVl)) . 
Indeed, for the main series, the effective charge Z\ dif- 
fers only slightly from the nominal charge of the cation 
(which indicates the predominantly ionic character of the 
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A-0 bond); for the other two groups, the charge is 
significantly larger, which indicates that the bonding be- 
comes more covalent in character.— 

The obtained results suggest that off-center impurity 
atoms can exist in solid solutions of titanates with the 
perovskite structure. Since the average interatomic dis- 
tance in these crystals is determined by the matrix, one 
can expect that the atoms for which ^xx{0, 0) has nega- 
tive or small positive values can be off-center. Therefore, 
it is possible that the ferroelectric phase transition in- 
duced by Ca, Cd, and Pb impurities in SrTiOa (Ref. [l^ 
is due to off-centering of these atoms. According to the 
EXAFS data, the Ba atoms in SrTiOs are on-center)^ 
but the Pb atoms in SrTiOa and BaTiOa can be off- 
center.'*'^ 

The results of this study differ somewhat from those 
obtained by Kvyatkovskiij^* according to which the mul- 
tiwell adiabatic potential occurred only for Mg and Zn 
atoms, whereas Cd atoms remained in the on-site posi- 
tions. The discrepancybetween our results is likely due 
to the fact that in Ref.|Tl|the calculations were performed 
for relatively small clusters in which the extended corre- 



lation of atomic motion cannot be correctly taken into ac- 
count (according to Refs. H and H, the correlation length 
can be as large as 20 A). 



VI. CONCLUSIONS 

The pseudopotentials constructed in this work have 
been used to calculate the phonon spectra of ^TiOs per- 
ovskite crystals within the density functional theory. All 
known results on the structural instability in these com- 
pounds were reproduced and the properties of new, pre- 
viously unknown systems were predicted. The analysis 
of the phonon spectra, the force constant matrix, and the 
eigenvectors of unstable TO phonons enabled to establish 
the relative contributions of the chain instability and off- 
centering of atoms to the appearance of ferroelectricity 
in these compounds. The main factors determining the 
possible off-centering of the A atoms are the geometric 
size of these atoms and the configuration of their outer 
electron shell. 
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